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Two body data alone cannot determine the potential uniquely, one needs three-body data 
as well. A method is presented here which simultaneously fits local or nonlocal potentials 
to two-body and three-body observablcs. The interaction of composite particles, due to 
the Pauli effect and the indistinguishability of the constituent particles, is genuinely 
nonlocal. As an example, we use a Pauli-correct nonlocal fish-bone type optical model 
for the a — a potential and derive the fitting parameters such that it reproduces the 
two-a and three-a experimental data. 
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1. Introduction 



The aim of inverse scattering methods is to determine the operator from the spec- 
trum. If we assume nonrelativistic dynamics, the inverse problem boils down to 
determining the potential. Mathematics alone, however, cannot determine the po- 
tential uniquely. A general potential operator, taken in coordinate representation, 
is a nonlocal potential 

(r\V\r') = V(r,r') . (1) 
So, if we write the Schrodingcr equation in the form 

n 2 ( d 2 i(i + iy 



2m V dr 2 



Vi{r) 



^(r) = Ei/>(r), (2) 



where m is the mass of the particle and I is the angular momentum, we assume 
implicitly that the potential is local 

V l {r,r') = V l {r)S{r-r'). (3) 



f 
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This is a serious approximation since, for example the interaction between compos- 
ite particles, due to the internal motion of the constituents, is a genuine nonlocal 
potential. 

The physical input to the inverse scattering procedure is the bound-state energy 
and the phase shift, or something equivalent. The bound state wave function is 
square integrable and 

ipir) ~ exp(— nr), (4) 

as r — ► oo, where k = y/2m/K 2 \e\ and e is the bound-state energy. If the potential 
is of short-range type the scattering wave function at large distances looks like 

ip(r) ~ exp(-ikr + lir/2 + 6), (5) 

where k = y/2m/h 2 E, E is the energy and S is the phase shift. So, the spectrum 
of the operator is contained in the asymptotic behavior of the wave function. Or, 
other way round, the spectrum determines only the asymptotic behavior of the 
wave function. The internal part is arbitrary, at least to some extent. 

To determine the potential unambiguously we need to know the internal part 
of the wave function as well. In order that we gain some information about the 
internal part of the wave function we may add a third particle to the system. The 
third particle feels not only the asymptotic part, but also the whole wave function. 
So, three-body observables are sensitive not only to the asymptotic part, but also to 
the whole two-body wave function. Therefore, to determine the potential uniquely, 
we need to study both two-body and three-body systems. 

The aim of this work is to present a quantum mechanical approximation method 
which can handle realistic local or nonlocal potentials with long range Coulomb tail 
in two- and three-body dynamics. For this purpose we developed an integral equa- 
tion approach. For two-body systems we solve the Lippmann-Schwinger integral 
equation, and for three-body systems we solve the Faddeev integral equations. Here 
we sketch the solution method, the details are given in preceding publication£H 

An argument against using nonlocal potentials is that it is too general and by 
using them we open a "Pandora's box". However, we often have some information 
about the internal structure of the composite particles which provides some frame- 
work for the nature of the nonlocality. Composite particles often made of identical 
elementary fermions. If two composite particles gets closer, the Pauli principle pre- 
vent the identical constituents from occupying the same quantum state. This is the 
main source of the nonlocality of the interaction at short distances. There are sev- 
eral cluster model inspired effective potentials which has nonlocal parts. Probably 
the most general and best founded model is the fish-bone model ^ This approach 
takes into account not only the complete Pauli blocking, but also it allows that 
some states are partly Pauli forbidden or Pauli suppressed. 

In this work we examine the a — a interaction in the fish-bone model. We 
determine the a — a potential such that it reproduces the S, D and G wave phase 
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shifts up to 20 MeV, the famous the E = 91 KeV, l = 0«-a resonant state, and 
the three-a binding energy. 

In Sec. II we present the Coulomb-Sturmian separable expansion approach to 
the two-body Lippmann-Schwinger equatiorW In Sec. Ill we extend the method for 
solving the Faddeev integral equations of the three-body problenP. The Faddeev 
equations with Coulomb potentials are rather complicated, so we restrict ourself 
to bound states and identical particles. Then, in Sec. IV, we reexamine the two-a 
and the three-a system and determine a new fish-bone-type a — a potential which 
provides a good description to the two-a and three-a observablcs. 



2. Two-body problem 

2.1. Coulomb-Sturmian separable expansion approach to 
Lippmann-Schwinger equation 

In this section we recapitulate the solution of the Lippmann Schwingcr equation 
by using the Coulomb-Sturmian separable expansion method. We suppose that the 
total Hamiltonian hi can be split into two terms 

h = hf + vi , (6) 

where vi is an asymptotically irrelevant short-range potential and hf is the Coulomb 
Hamiltonian 



■ c _ h\ 2 f d 2 + , Z 



2m \ dr 2 



r 



(7) 



where m is the reduced mass, I is the angular momentum and Z is the strength of 
the Coulomb potential. 

We denote the Coulomb Green's operator by gf"(z) = (z—hf)^ 1 and the Green's 
operator of the total Hamiltonian by gi{z) = (z — hi) -1 . They are related via the 
resolvent relation 



gi{z) =gf{z)+g\ ?{z)v igi {z) . (8) 

(- 
i 



The scattering wave function satisfies the inhomogeneous Lippmann 



Schwinger equation 

\ip\ +) ) = W?) + g?{E + \0)vM\ +) ) , (9) 

while the bound- and resonant-state wave function satisfies the homogeneous 
Lippmann-Schwinger equation 

\^ l )=g?(E)v l \^ l ) (10) 

at negative real and complex E energies, respectively. 

We solve these equations in a unified way by approximating only the short range 
potential vi. We use Coulomb-Sturmian (CS) functions 

(r\nl;b) = ( - ^ + ) " exp(-fer)(26r) i+1 4 2 ' +1 )(26r) , (11) 
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where L denotes the Laguerre polynomial. The CS functions form a bi-orthonormal 
basis. With (r\nl; b) = (r\nl; b)/r we have 

(nfl{b\nl, b) = (ril; b\rdfb) = 5 nn ,. (12) 

To approximate the potential operator we use a " skew" form of the completeness 
relation 1; — limjv-^oo if, where 

N 

lf= E HibgKQ-^nmiml;^] (13) 

n,m— 

with (0) mn — ((ml; b v \nl\b g )). It is obvious that this double sum is also a possible 
expression for unity. If b v = b g , Eq. ( 13 ) falls back to the usual completeness relation. 
By adopting different values for b v and b g one can achieve a faster convergence in 
the separable approximation^. 

In Eqs. ([9| and ( |To| ) the term vi\4>i) is square integrable, and belongs to the 
Hilbcrt space. Therefore, it can be approximated by 

N 

w,M) « Ifv^) « l? Vl lf |^) = «/>,) = E \ nl ^ 9 ) V (14) 

n.n' 

with finite AT, where 

N 

V= E (0 _1 )„ m (TO/;6 1 ,|fi|m'?;6 t ,)(0~ 1 ) m /„/ . (15) 

m,m'— 

In general, the matrix elements (mZ; &t,|i>z|m'Z; &„) have to be calculated numeri- 
cally. Note that the potential is sandwiched between CS states with parameter 
b v and the potential operator becomes a linear combination of CS ket-bra opera- 
tors \nl;bg}(n'l;b g \. This approximation is called separable expansion because the 
potential operator after approximation appears in the form 

N 

(rK\r') = J2 ( r K & 9 > !W WsV) , (16) 

n,n'=0 

i.e. the dependence on r and r' is separated. 



With this separable potential Eqs. (J9J) and (10 1 become 

N 



= \<P?) + E 9?(E + i0)\nl;b g ) v nn , (n'l;b g \4 +) ) , (17) 

n,n'=0 

and 

JV 

IV*) = E 9?(E)\nl;b g ) v nn , (n'l;b g \^) , (18) 



respectively. To derive equations for the unknown coefficients = (n'l; b g \ip[ + ^) 
and tj) — (n'l;b g \ipi), we have to act with states (n"l;b g \ from the left. Then, the 
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following inhomogeneous and homogeneous algebraic equations are obtained for 
scattering and bound-state problems, respectively: 

[(gf(E + iO))" 1 - = (gf(E + iQ))"Vf , (19) 

and 

[(gf(E))-'-v i }^ = 0. (20) 

Here we have g^ = (nl; b g \gf \ (n'l; b g ) , the matrix elements of the Coulomb Green's 

operator, and the overlap of the CS and Coulomb functions <f>? = (n'l; b g \4>f). The 
reason that we have adopted the CS basis is that, in that basis, these quantities 
can be calculated analyticalljEE The homogeneous equation (20 1 is solvable if and 
only if 

det[(gf(£;))- 1 -^ = (21) 
holds, which is an implicit nonlinear equation for the bound- and resonant-state 



energies. As far as the scattering states arc concerned, the solution of (19) pro- 
vides the overlap (nl; b g \^pi). From this quantity any scattering information can be 
inferred, for example the Coulomb-modified scattering amplitude reads 

a l = (^- ) \v l \4 +) )=^v l ^ l +) , (22) 
which is related to the Coulomb-modified short-range phase shift 5i through 

ai = ^exp(i(2?7 ; +(5 ; ))sin<5 / , (23) 

where rji = argT(Z + 17 + 1) is the Coulomb phase shift. 

In this approach, the potential enters into the procedure through its CS matrix 
elements. So, the method works for any potential, local or nonlocal, energy depen- 
dent, complex, etc., as long as we can evaluate its CS matrix elements somehow. 
Numerical examples show that the method is also very efficient 1 . Having the opti- 
mal values for b v and b g , one need about N = 15 — 20 basis states to achieve 6 — 7 
digits accuracy for the bound state energies and for the phase shifts over the entire 
spectrum. 



3. Three-body problem 
3.1. The Faddeev equations 

Here we consider three identical particles interacting with repulsive Coulomb-like 
potentials and we restrict ourself to bound states. The method is applicable also 
for attractive Coulomb potentials and for scattering and resonant states as welP. 
The Hamiltonian is given by 



H = H° + Vl + v 2 + v 3l 



(24) 
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where H° is the three-body kinetic energy operator and Vi denotes the long-range 
Coulomb-like potential of each subsystem i = 1,2,3. We use the usual configuration- 
space Jacobi coordinates: X\ is the distance between the pair (2, 3) and yi is the 
distance between the center of mass of the pair (2, 3) and the particle 1. 

A Coulomb potential modifies the character of the asymptotic motion, there- 
fore it should be treated very much like the kinetic energy operator. We split the 
potential into two parts, a short-range and a Coulomb part 

v i = vV+v?. (25) 

We define the long-range Hamiltonian by 

H d) = h° + + v$ + v$, 

and the three-body Hamiltonian takes the form 

H=H^+v{ s) +vi s) + vi s) . (27) 

This Hamiltonian looks like an ordinary three-body Hamiltonian with short-range 
interactions. We solve the Schrodinger equation 

H\V)=E\$) (28) 

by using the Faddeev method. We split the wave function into three components 

|*) = + \ifa) + \ips), (29) 

and the components should satisfy the set of differential equations 



(E-Hl l) )\i> 1 )=vl s) (\i> 2 ) 


+ W) 


(30a) 


(E-Hi l) )\^) = vi s \\^) 


+ |^3» 


(30b) 


(£-i/f)i^) = 4 s) (hM 


+ M*»> 


(30c) 


where Hf — + . By adding these three equations we 


get back the original 


Schrodinger equation. 






By inverting the left hand side of (|30| we get a set of intef 


;ral eqautions 


\1>i) = G«\E)vi'\\ih)-{ 


-1^3)) 


(31a) 


|^ 2 > = Gf (£)4 s) (|V>i>4 


" |^3» 


(31b) 


\i> 3 ) = G« ) (E)vi s) (\i> 1 )^ 




(31c) 


where G ( - ] {E) = {E-H^y 1 . 







If particles 1, 2 and 3 are identical, then ipi, in its natural Jacobi coordinate sys- 
tem {x\, yi}, looks like ip2 in its natural Jacobi coordinate system {^2, 2/2} and tps in 
its natural Jacobi coordinate system {£3,2/3}. On the other hand, by interchanging 
particles 2 and 3 we have 



(26) 



^23 l^l) =P\lpl), 



(32) 
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where p = 1 for bosons and p = — 1 for fermions. Building this information into the 
formalism we arrive at a single integral equation 

|^} = 2GfM S) ?Wi), (33) 

where V123 = V12P23 is the operator for cyclic permutation of all three particles 
^'123 IV^i) — l^)- We should notice that so far no approximation has been made, 
and even though this integral equation has only one component, yet it gives a full 
account for the asymptotic and symmetry properties of the system. 



3.2. Coulomb- Sturmian expansion 

Since the three-body Hilbert space is a direct product of two-body Hilbert spaces, 
an appropriate basis is the bipolar basis, which can be defined as an angular- 
momentum-coupled direct product of the two-body bases, 

\nvl\;b x by) a = \nl;b x ) a ig)\i>\;by) a , (n,u = 0,1,2, .. .), (34) 

where \nl; b x ) a and \v\;b v ) a are associated with the coordinates x a and y Q , re- 
spectively. With this basis the completeness relation takes the form (with angular 
momentum summation implicitly included) 

jv 

1 lim y \nvl\;b x b y ) a a {nvl\;b x by\ = lim 1^, (35) 

AT— >oo ^ — 4 TV— >00 

where (xy\nvl\\b x b v ) — (xy\nisl\;b x b y )/(xy). 

Similarly to the two-body case, v^V^l^i) is square integrable, therefore we 
can approximate 



v{ s) V 123 \^) = lim If v[ s) V 123 l» » If v{^P 123 l? |Vi> 

N^oc 

N _ _ 

J2 Invl^b^tV^in'v'l'X'-^by]^), (36) 



n.v.n' v'— 



where V_ 1 = y^{nl; b x \n'l; b x ) i(nvl\; b x b y \n' V 'V A' '; b x b y )2- The completeness of the 
CS basis guarantees the convergence of the expansion with increasing N and angular 
momentum channels. 

Now, by applying the bra (n"v"l"\"\ b x b y \ on Eq. (33 1 from the left, the solution 
of the homogeneous Faddeev equation turns into the solution of a matrix equation 
for the component vector ip = \{nvlX; b x b y \ip\) 

± x = 2GfVi^, (37) 

where 

Gf = 1 (niJW x by\GP\n>vn ; \Xb x b y ) v (38) 
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This homogeneous algebraic equation is solvable if only if the determinant is zero: 

det[(GiV X - 2Zi] = 0. (39) 

The operator Crf 1 is the resolvent of a complicated three-body Coulomb Hamil- 
tonian. This Hamiltonian, however can support only one kind of asymptotic chan- 
nel, when particles 2 and 3 are close and particle 1 is at infinity. The asymptotic 
Hamiltonian of is 

H 1 = H° + v 1 +u 1 , (40) 

where is the channel Coulomb potential. This potential is the asymptotic part 
of v§ + as particle 1 is separated from the pair (2, 3). It looks like the charges 
of particles 2 and 3 are concentrated in their center of mass 



c 



Z 1 (Z 2 + Z 3 )/y 1 , (41) 



where Zi is the charge of the particles. 

For Hamiltonians which have one kind of asymptotic Hamiltonian, a single 
Lippmann-Schwinger equation provides a unique solution. Thus 

Gf\z) = G t (z) + G 1 {z)U^G ( (\z), (42) 

where G\(z) = {z — H\)~ l and — v% + — u^. To solve the Faddeev equation, 
we need G^ between finite number of square integrable CS basis states. We make 
a separable approximation on U\, and with the help of matrix elements U_ 1 — 
i{nvl\;b x b y \Ui\n' v'V \' ;b x by)i, we get 

(G^)- 1 ^)- 1 -^, (43) 

where G ± = i{nvl\;b x b v \Gi\n'i>'l'\';b x by)i. This later matrix element of U\ can 
always be evaluated numerically. 

To calculate the matrix elements G 1 we utilize the Dunford- Taylor functional 
calculus. If h is a selfadjoint operator, then an analytic function / of h is given by 
a contour integral 

f(h) = ^~ <f dzf(z)(z-h)-\ (44) 

where the contour C goes around the spectrum of h such that / is analytic on the 
area encircled by C. 

The three-particle free Hamiltonian can be written as a sum of two-particle free 
Hamiltonians 

H° = h Xi + h° yi . (45) 



Thus the Hamiltonian Hi of Eq. ( 40 1 appears as a sum of two two-body Hamilto- 



nians acting on different coordinates 

Hi=h Xl +h yi) (46) 
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with h Xl = h xi + Vi(x±) and h Vl = hy x + (yi), which, of course, commute. The 
Green's operator 

G 1 {E) = {E-h Xl -h yi )- 1 (47) 
is a function of the selfadjoint operator h yi . So, 

Gi(E) = -L I dz (E-h^-z^iz-hhJ- 1 = ^.l dzg Xl (E-z)g yi (z),(48) 
^ 7ri Jc ^ 7Tl Jc 

i.e. we can calculate the three-body resolvent as a convolution integral of two-body 
resolvent operators. 

In this work we concentrate on bound states of the three-a system. Therefore 
E < 0, h yi is a repulsive Coulomb Hamiltonian, and although h Xl contains a 
short-range potential, it does not support bound states. Figure I shows the analytic 



structure on the integrand ( 48 ) and the contour which encircles the spectrum of h yi 
without penetrating into the spectrum of h Xl . In Fig. II the contour C is deformed 
analytically. This is a more advantageous contour for numerical calculations since 
the matrix elements of Green's operators falls off quickly and smoothly in this 
complex imaginary direction. 



Fig. 1. A schematic illustration of the analytic structure of the integrand ( |48| . If E < 0, the 
two branch-cuts are well separated. The contour C goes around the spectrum of h yi without 
penetrating into the spectrum of h xi . 



g 



C 



Fig. 2. The same as in Fig. 1. Here the contour is deformed analytically to achieve a smother 
integrand. 
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Table 1. L = three-a binding 
energy as a function of subsys- 
tem angular momentum l m 

ax in 

case if Ali-Bodmer (AB) poten- 
tial, fish-bone potential of Kirchcr 
and Schmid (FB-1) and the re- 
sults of this work (FB-2). 



Imax 


AB 


FB-1 


FB-2 


2 


-1.70 


0.057 


-0.313 


4 


-2.26 


-15.47 


-7.112 


6 


-2.28 


-15.63 


-7.273 


8 


-2.28 


-15.63 


-7.275 



4. The a. — a. potential 

Two- and three-a systems have been subject to a very intensive study over the past 
decades. We just refer here to some recent investigations ranging from local plus 
three-body potential models^, through orthogonality condition modeP, to genuine 
resonating group calculations^. 

4.1. The Ali-Bodmer potential 

The "standard model" for the a — a interaction is the Ali-Bodmer potentiaP. There 
are various parameterizations but, in general, in this model, the potential is local 
and partial- wave dependent. We use the parametrization 

{-150cxp(-0.5r 2 ) + 1050exp(-0.8r 2 ) if I = 
-150exp(-0.5r 2 ) + 640exp(-0.8r 2 ) if I = 2 (49) 
-150exp(-0.5r 2 ) if I > 4 even, 

which is shown in Fig. III. We can see that the potential has a very strong angular 
momentum dependence. 

Table I shows the binding energy of the three-a system with total angular 
momentum L = as a function of maximal subsystems angular momenta l max — 
\nax- We used CS function up to N — 25, which provided us with about five digits 
accuracy. If we take the I = A = angular momentum channel only, we do not get 
binding. We can also sec that angular momentum channels I = A > 4 have very 
little contributions to the binding energy. The experimental binding energy of the 
three-a (C12) system is E 3a = —7.275 MeV, so the Ali-Bodmer potential does not 
provide enough binding. 

4.2. The fish-bone optical model 

The fish-bone model is motivated by the cluster model. In the resonating group 
model the total wave function is an antisymmetrized product of the cluster $ and 
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\ 
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\ 




Fig. 3. The Ali-Bormer a — a potential for I = 0, I = 2 and I = 4. 



the inter-cluster \ relative states 

I*) = |M*x}>. 



(50) 



The state <i>, which is supposed to be known in this model, describes the inter- 
nal properties of the clusters, including spin and isospin structure. The unknown 
relative motion state x is determined from the variational anzatz 



{®5 X \A(H - E)A\$x) = Q. (51) 

This anzatz results in a rather complicated equation for which were possible 
to solve only by using serious approximations on $ and on the interaction of the 
particles. Based on this resonating group cluster model, several models has been 
proposed to describe the motion and interaction of composite particles. Probably 
the most elaborated one is the fish-bone optical model proposed by SchmicP. 
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In the fish-bone model the interaction of the two-body subsystem is given by 



V; = vi - ) \u Li }(ui^\(h° + vi- e/,i)|«i, 3 -)Afi^(«j j-|, 



(52) 



where Z refers to partial wave, h is the kinetic energy and vi is a local potential. 
The states |uz,i) are eigenstates of the norm operator, 

($f*|./4|$u,,<} = (l-r? M )(r1u M ), (53) 

where r is the center of mass distance of the two clusters. If the relative motion is 
forbidden by Pauli principle then (^r^l^u;,,) = 0, and r\i t i — 1. The r\\^ eigenvalues 
are ordered such that \t)h\ > The matrix M is then given by 

. . (54) 

[(1-^.)(1_^)]V 

where 77;^ = if rjy = 1 and rjij — r]i t i otherwise. Or, in matrix form, if we have 
one Pauli forbidden state, 

/I 1 



1 



1-7- r— — ^1/2, rf*>Jj 



A/, 



1 

1 1 - 
1 1 - 







l— m, 
i-»n. 



1 


-»7(,2 




-171,3 


'1 


-»)1,2 




->7l,3 



1 


-»)!,2 




-»7i,4 


1 


-»7l,2 




-»)i,4 



1 - 



l-»?i,2 
1— m,4 







(55) 



V: •••/ 

which exhibits a fish-bone-like structure, where the name of the model comes from. 
In this model the Pauli-forbidden states become eigenstates at e energy. By choos- 
ing e as large positive, they become bound states at large positive energy, and 
thus disappear from the physically relevant part of the spectrum. There are sev- 
eral versions of the fish-bone model which differ in off-shell transformations, i.e. in 
transformations which effect the internal part of the wave function and leave the 
asymptotic part, and the spectrum, unchanged. This version of the model minimizes 
the three-body potential, which is therefore neglected. 

We assume that in the a particles the nucleons are in 0s states in an os- 
cillator well of width parameter a. Then the norm kernel eigenvalues are also 
harmonic oscillator functions with the same width parameter and the eigenval- 
ues are knowrP! r] 0>i = 1,1,1/4,1/16,1/64,..., 773,, = 1,1/4,1/16,1/64,... and 
774 j = 1/4, 1/16,1/64, . . .. So, in the I — relative motion channel there are two 
Pauli-forbidden states, in / = 2 there is one, an in I — 4 and higher channels there 
are none. The decreasing value of rj indicates that those harmonic oscillator sates 
in the relative motion are less and less suppressed by the Pauli principle. For the 
e parameter of the fish-bone model, which aim is to remove the Pauli-forbidden 
states, we took e = 200000 MeV. In this range of e, the dependence of the results 
was beyond the fifth significant digit. 
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In an earlier work by Kircher and SchmicP^ a fish-bone-type potential was deter- 
mined from two-a scattering data. The width parameter was fixed to a = 0.55fm -2 
and the the local potential was given by 

4e 2 / [2a \ 

vi(r) = v exp(-/3r 2 ) H erf \ —r , (56) 

r \\ 3 ; 

where v = -108.41998MeV and = 0.18898fm~ 2 . While this potential provides a 
reasonably good fit to I = 0, I = 2 and I = 4 partial wave phase shifts, it seriously 
overbinds the three-a system (see FB-1 in Table I). One may conclude that there is 
a need for three-body potential. This was the choice Oryu and Kamada^ adopted. 
They added a phenomenological three-body potential to the fish-bone potential 
of Kircher and Schmid and found that a huge three-body potential is needed to 
reproduce the experimental data. But, our Faddeev calculations reveal that the 
I = 4 partial wave is very important to the three-a binding and, for this partial 
wave, the fit to experimental data is not so stellar. So we concluded, that it may be 
possible to improve the agreement in the I — 4 partial wave and achieve a better 
description for the three-a binding energy. 

As a local potential we took two Gaussians plus screened Coulomb potential 

vi(r) = vi exp(-/3ir 2 ) + v 2 exp(-/3 2 r 2 ) + ^-erf ^\/y r ^j > ( 57 ) 

where V\, fli, v 2 , $i and a are fitting parameters. In the fitting procedure we in- 
corporated the famous s Be, I — resonance state at E^ p — (0.0916 — 0.000003i) 
MeV, the 12 C three-a ground state energy E^ p = -7.275 MeV, and the I = 0, 
I = 2 and I = 4 low energy phase shifts. With parameters vi = —120.30683493 
MeV, Pi = 0.20206127 fm~ 2 , v 2 = 49.06187648 MeV, (3 2 = 0.76601097 fm" 2 and 
a = 0.64874009 fm~ 2 we achieved a perfect fit. For the I — two-body resonance 
state we get E 2b = 0.09161 - 0.00000303i MeV, and for the three-body ground 
state Est, = —7.27502 MeV. The fit to the phase shifts is shown on Fig 4. We can 
see that the fit could hardly be any better. Fig 5 shows the local part of the fish- 
bone potential. Notice that unlike with the Ali-Bodmer potential, we achieved this 
agreement by using the same potential in all partial waves. 

Having this new a — a fish-bone potential from the fitting procedure, we also 
calculated the first excited state of the three-a system. This state is a resonant 
state, and we got E^ s — (0.54 — 0.0005i) MeV, which is again very close to the 
experimental value. 



5. Summary 

In order that we can incorporate the Pauli effect we have to leave the local potential 
models and have to use nonlocal potentials. We solved the underlying Lippmann- 
Schwinger and Faddeev integral equations by performing a separable expansion 
of the potential in Coulomb-Sturmian basis. This expansion scheme can handle 
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Center of mass energy (in MeV) 



Fig. 4. Fit to the experimental I = 0, I = 2 and I = 4 phase shifts. 



practically any kind of Coulomb plus nonlocal short-range potentials. The Faddeev 
equations have a further advantage that they allow a very detailed channel-by- 
channel analysis of the interaction. The method is also very efficient numerically. To 
solve the Faddeev equation for five angular momentum channels takes a few minutes 
on our Mac PC, and to fit the fish-bone potential parameters took a couple of hours. 
We believe that we possess the right model and appropriate tools for investigating 
the interaction of composite particles. 

It is truly remarkable that we got a unified description of the two-a I = 
resonant state and the I = 0, I = 2 and I = 4 phase shifts, as well as the threc-a 
binding energy. This was achieved without resorting to three-body potentials and 
adding explicit angular momentum dependence to the local potential. The results 
show that the angular momentum dependence comes form the Pauli effect. If the 
Pauli effect is properly incorporated, like in the fish-bone model, and the fit is done 
in such a way that both two- and three-body data are incorporated, there is no need 
for additional angular momentum dependence and there is no room for a strong 
three-body potential. We can see that the proper inclusion of the Pauli principle 
simplifies the model potential of composite particles. 
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Fig. 5. The local part of the fish-bone potential uj. 
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